		twoway	(scatter $p vars1 if dataset == "bs_decomposition_SS_ProdR_new_kl", mcolor($cl) msymbol($sym)) ///
				(rcap lows_$p highs_$p vars1 if id == 2 & dataset == "bs_decomposition_SS_ProdR_new_kl", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars1 if id == 1 & dataset == "bs_decomposition_SS_ProdR_new_kl", lwidth(thin) color(gs6) ) ///
				(scatter $p vars2 if dataset == "bs_decomposition_SS_ProdR_new_kl_3lags", mcolor($cl) msymbol(Oh)) ///
				(rcap lows_$p highs_$p vars2 if id == 2 & dataset == "bs_decomposition_SS_ProdR_new_kl_3lags", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars2 if id == 1 & dataset == "bs_decomposition_SS_ProdR_new_kl_3lags", lwidth(thin) color(gs6) ) ///
				(scatter $p vars3 if dataset == "bs_decomposition_SS_ProdQ_new_kl", mcolor($cl) msymbol(Dh)) ///
				(rcap lows_$p highs_$p vars3 if id == 2 & dataset == "bs_decomposition_SS_ProdQ_new_kl", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars3 if id == 1 & dataset == "bs_decomposition_SS_ProdQ_new_kl", lwidth(thin) color(gs6) ) ///
				(scatter $p vars4 if dataset == "bs_decomposition_SS_CS_kl_nic4", mcolor($cl) msymbol(Th)) ///
				(rcap lows_$p highs_$p vars4 if id == 2 & dataset == "bs_decomposition_SS_CS_kl_nic4", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars4 if id == 1 & dataset == "bs_decomposition_SS_CS_kl_nic4", lwidth(thin) color(gs6) ) ///
				(scatter $p vars5 if dataset == "bs_decomposition_SS_CS_kl_nic2_by_year", mcolor($cl) msymbol(Sh)) ///
				(rcap lows_$p highs_$p vars5 if id == 2 & dataset == "bs_decomposition_SS_CS_kl_nic2_by_year", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars5 if id == 1 & dataset == "bs_decomposition_SS_CS_kl_nic2_by_year", lwidth(thin) color(gs6) ) ///
				, xlabel(2 "R. TFP" 4 "M. Demand" 6 "M. TFPQ" 8 "Covariance" 10 "Scope" 12 "Adding" 14 "Dropping", angle(35)) ///
				ylabel(,grid) ///
				legend(order(1 "KL" 4 "KL 3Lags" 7 "Q KL" 10 "CS KL NIC4" 13 "CS KL NIC2 by Year") size(medlarge) rows(2) region(lcolor(white)) position(6) ring(1)) scheme(s1color) 
				*title("GO $p") 
	

			graph export "$figures/robustness_SVA_$p.png", as(png) replace width(1800) height(1000)
